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Abstract 

The light pseudoscalar mesons play a twofold role: they may or have to be regarded both as 
low-lying bound states of the fundamental degrees of freedom of quantum chromodynamics 
as well as the (pseudo-) Goldstone bosons of the spontaneously broken chiral symmetries of 
quantum chromodynamics. We interrelate these aspects in a single quantum-held-theoretic 
approach relying on the Bethe-Salpeter formalism in instantaneous approximation by very 
simple means: the shape of the pseudoscalar-meson Bethe-Salpeter wave function dictated 
by chiral symmetry is used in Bethe-Salpeter equations for bound states of vanishing mass, 
in order to deduce analytically the interactions which govern the bound states under study. 
In this way, we obtain exact Bethe-Salpeter solutions for pseudoscalar mesons, in the sense 
of establishing the rigorous relationship between, on the one hand, the relevant interactions 
and, on the other hand, the Bethe-Salpeter amplitudes that characterize the bound states. 
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1 Introduction 

Light pseudoscalar mesons may be understood as bound states of a quark and an antiquark 
but they have to be interpreted also as (almost) massless (pseudo-) Goldstone bosons of the 
spontaneously broken chiral symmetries of quantum chromodynamics, the theory of strong 
interactions, whence their description still poses a challenge to theoretical particle physics. 
Within quantum field theory, the Bethe-Salpeter framework provides a PoincarGcovariant 
approach to bound states [1]. Some of its inherent obstacles can be avoided by restriction to 
three-dimensional, e.g. instantaneous, reductions, such as Salpeter’s equation [2], Recently, 
by inversion of this bound-state problem, we showed that, under favourable circumstances, 
the underlying interaction potential can be retrieved from the Salpeter solutions [3]. In this 
paper, we apply the inversion technique of Ref. [3] ^ to pseudoscalar Goldstone-type mesons 
(like the pion), to see how the strong interactions enter such kind of bound-state equation.^ 
The outline of this paper is as follows. In Sec. 2, we recall briefly and thus in a somewhat 
symbolic notation only those aspects of the Bethe-Salpeter formalism that are of relevance 
for the subsequent discussion. In Sec. 3, we sketch the, in fact, not particularly complicated 
concepts behind our inversion procedures. In Sec. 4, in an attempt to get the most out of it, 
we squeeze dry the relevant results emerging from Euclidean-space-based analyses utilizing 
the Dyson-Schwinger framework, in order to extract from these well-motivated conjectures 
how the bound-state amplitudes forming the starting point of the inversion might look like. 
Armed with these insights, it is just one small step for (a) man to recover, in Sec. 5, for both 
zero (Subsec. 5.1) and non-zero (Subsec. 5.2) equal mass of the two (anti-) quarks bound to 
Goldstone bosons the basic interactions. Finally, Sec. 6 is devoted to summarizing remarks. 


2 Bethe-Salpeter Formalism in Instantaneous Limit 

An inevitable prerequisite of the application of our inversion approach developed in Ref. [3] 
to quark-antiquark bound states of Goldstone type is clearly the sufficient simplification of 
any quantum-field-theoretic description of bound states. Let us thus recall its crucial steps. 

The Bethe-Salpeter formalism [1] describes the features of bound states of fundamental 
degrees of freedom of one’s quantum held theory in use by the solutions of the homogeneous 
Bethe-Salpeter equation, deduced in form of pairs (M, <h) of its (discrete) bound-state mass 
eigenvalues M and associated eigenstates, represented by Bethe-Salpeter amplitudes <h. In 
momentum-space representation, each such Bethe-Salpeter amplitude <h(p, P) encodes the 
distribution of the relative momenta p of the two constituents of the respective bound state 
of total momentum P. The Bethe-Salpeter equation relates this Bethe-Salpeter amplitude 
<h(p, P), for bound-state constituents having momenta p*, i = 1,2, to their full propagators 
Si{pi) and an integral kernel K{p, q, P) that encompasses all interactions of these particles: 

Hp,P) = / dUK{p,q,P)^q,P)S2{-p2) • ( 1 ) 

^In Ref. [3], we focused (mainly) to the construction of exact solutions to the reduced Salpeter equation. 
The reduced Salpeter equation emerges from the (full) Salpeter equation when assuming the systems under 
consideration to be composed of weakly bound semirelativistic heavy constituents. In contrast, the present 
analysis deals with the full Salpeter equation, which does not need any such requirement. Hence, the rather 
severe constraints on the nature of the bound states do not apply in the case of the (full) Salpeter equation. 

^ We would like to convey in this way our gratitude to the referee of Ref. [3] for instigating these analyses 
by expressing great interest in the application of the formalism elaborated in Ref. [3] to the case of the pion. 
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In order to render this Bethe-Salpeter formalism accessible to our inversion techniques, 
we consider some so-called three-dimensional reduction of this framework. The by far most 
popular among such class of approximations is the static limit, encountered if assuming the 
bound-state constituents, in their center-of-momentum frame dehned by p = pi = —p 2 , to 
interact instantaneously. The kernel then depends only on the relative three-momenta p, q: 


K{p,q,P)=K{p,q) . 


If the propagators Si{pi) do not involve any non-trivial dependence on the time component 
Pq of the relative momentum p, the Bethe-Salpeter equation (1) becomes, upon integration 
over Pq, a generalized [4] instantaneous Bethe-Salpeter equation for the Salpeter amplitude 

0(P) = ^/dPo^(p) • (2) 


An example of such kind of equation has been derived in Refs. [4,5] and explored in Ref. [6]. 

The easiest way to accomplish the desired triviality of the po-dependence of any fermion 
propagator is to follow Salpeter [2] by assuming each bound-state constituent to propagate 
freely and, accordingly, approximating in Eq. (1) each full propagator S'j(p) by its free form 
Sifi{p, rrii), with, however, an effective mass m* subsuming all dynamical self-energy effects: 


Si,o{p, rrii) 


i 

f — rrii+ \e 


f + rrii 

9 9 i ^ ’ 

p^ — mf + le 


^ = ^10, f = l,2. 


Expressed in terms of the one-particle energy Ei{p), one-particle Dirac Hamiltonian Hi{p), 
and energy projection operators Ef{p) for positive and negative energy of particle i = 1,2, 

Ei{p) = sjp‘^ + ml , Hi{p) = 7o (7 ■ p + mi) , Af{p) = ’ 

the result of these simplihcations of the Bethe-Salpeter equation (1) is (upon application of 
contour integration in the complex-po plane and residue theorem) the Salpeter equation [2] 


A( f M^(Pi)7o [A:(p,q)0(g)]7oA2 (P 2 ) 

J {2ny\ Pq - E,{p,) - EM 

jPi) 7o [K{p, q) 0(g)] 7o (P 2 ) 
Pq + Ei{pi) p E2{P2) 


( 3 ) 


Its specihc projector structure subjects the Salpeter amplitude (2) to the crucial constraint 


A+ (pi) 0(p) A^(P 2 ) = Ai (pi) 0(p) A 2 (P 2 ) = 0 . (4) 


Counting the number of basis elements of the Dirac algebra, the Salpeter amplitude can 
have, in principle, at most 16 independent components. The constraints (4), however, halve 
this number: the most general Salpeter amplitude 4>{p) has eight independent components. 
The Salpeter amplitude of bound states of a spin-i fermion and a spin-| antifermion whose 
spin quantum numbers add up to zero has merely two independent components, henceforth 
labelled (pi{p) and 9 ^ 2 (p); pseudoscalar bound states are just that special case of this where 
the relative orbital angular momentum i of these bound-state constituents vanishes as well. 
The constraint (4) enforces, as general form of any such [CP = —1) Salpeter amplitude [7], 


0(P) 




75 • 


( 5 ) 
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The interaction kernel K{p, q) in Salpeter’s equation can be written in form of a sum of 
products of tensor products Ti^T 2 of Dirac matrices T 12 , defining the Lorentz structure of 
the effective couplings of the fermions, and associated Lorentz-scalar potentials generically 
labelled D(p, q). Here, we find reasonable to assume the equality ri = r 2 = rofri and r 2 : 

[K{p,q)(t){q)\ = ^VT{p,q)T (t){q)V . ( 6 ) 

r 

Convolution nature and spherical symmetry of all functions Vr(p, q) = Hr((p—— and 
hence of the entire kernel K{p, q) = K{{p—q)‘^) — imply that the Fourier transform of any 
such function Vr{p, q) is a conhguration-space central potential Vr(r), r = |a:;|. Splitting off 
all dependence on the angular variables then converts Salpeter’s equation (3) to a system of 
coupled equations for the radial factors of the independent components [7]. Suppressing the 
index F, any potential V (r) enters in these radial equations by its Fourier-Bessel transform 

00 

VL{p,q) = 87rJdrr‘^jL{pr)jL{qr)V{r), p = \p\ , Q = \q\ , L = 0 , 1 , 2 ,..., 

0 

given in terms of the spherical Bessel functions of the hrst kind [ 8 ] jniz), n = 0, ±1, ±2,.... 

Let us now zoom in onto our actual targets, i.e., the light pseudoscalar mesons: ordinary 
— in contrast to exotic^ — mesons that may be understood as bound states of a light quark 
and a light antiquark. Ignoring flavour violation enforces equality of the light-quark masses 
mi = m 2 = m and hence of the free energies: Ei{p) = 7^2 (p) = E{p) = ^/p'^ + m^, p = \p\. 

In order to dehne unambiguously the instantaneous Bethe-Salpeter equation we intend 
to invert, we need to specify the Lorentz structure of all tensor products F( 8 )F of generalized 
Dirac matrices F in the interaction kernel ( 6 ). Following Refs. [9] (which form examples of a 
phenomenologically acceptable analysis of the quark-antiquark bound state spectrum), we 
choose F 0 F to be the unique (whence the right-hand side of Eq. ( 6 ) reduces to only a single 
term) sophisticated linear combination of scalar, pseudoscalar, and vector Dirac structures 

F 0 F = ^ ( 7 ^ 0 7 ^ -F 75 0 75 - 1 0 1) , (7) 

which has the distinctive feature of Fierz symmetry, i.e., invariance under rearrangement of 
Dirac helds ffkix), k = 1,..., 4: ffii^x) T ■^ 2 ( 2 ^) F 'ip 4 {x) = 'ipi{x) F 'ip 4 {x) F 'ip 2 {x). 
With the educated guess (7) for the kernel, the Salpeter equation (3) for spin-singlet bound 
states becomes equivalent to a set of coupled eigenvalue equations for the two radial factors 
ipi{p) and ip 2 {p) of the independent components of the Salpeter amplitude (5) [10, Sec. IX]: 

00 

2E{p)(p2{p)+ 2 j ^^Vo{p,q)^2{q) = M ^i{p) , 

0 

2E{p)^i{p) = M(p2{p) ■ 

^For a fermion-antifermion bound state, its parity P and — if this bound state is composed of a fermion 
and its associated antiparticle and therefore exhibits a well-defined behaviour under charge conjugation — 
its charge-conjugation parity C are related to the relative orbital angular momentum £ and the total spin S 
of the bound-state constituents according to P = (—1)^+^ and C = (—1)^+'^; for any such bound state with 
total spin J conceivable quantum-number assignments are = 0++, 0 1^“'', , 1 , 2++,.... 

Mesons carrying any such assignment may be quark-antiquark bound states; their complement, those with 
an assignment , 0 ,1 2“' ,3 4^ ,..., must be non-qq and belong to the exotic mesons. 
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The eigenvalues of this eigenvalue problem are the possible bound-state masses M = 

The hrst of these coupled relations is an integral equation encompassing all the information 
about the interactions, whereas the second one is of merely algebraic nature. For the case of 
interest here, i.e., for vanishing (Goldstone-boson) mass M = 0, the two relations decouple: 

• The second relation implies that one of the Salpeter components vanishes: 9 ?i(p) = 0. 

• The only non-zero Salpeter component, 9 ^ 2 (p), must satisfy the bound-state equation 

00 

E{p) P 2 {.p) + j Vo{p, q) p 2 {.q) = 0 , (8) 

0 

which, because of particularly fortunate circumstances, is equivalent to what is called 
the spinless Salpeter equation. (For reviews on the latter, consult, e.g.. Refs. [11-15].) 


3 Inversion Procedure 

The actual goal of our inversion technique [3] is, for a given bound-state Salpeter amplitude 
^(p), the extraction of the underlying conhguration-space potential V (r) from the relevant 
bound-state equation carved out within a sufficiently simplihed Bethe-Salpeter formalism. 
Preferably, this extraction should be accomplished (as far as possible) by analytical means. 

As has been shown in Ref. [3] by various examples, this main goal can be easily achieved 
if representing, by application of a Fourier transformation, the bound-state equation in use 
in conhguration space. The Fourier transformation in three dimensions of any purely radial 
function reduces to its (L = 0) Fourier-Bessel transformation involving the n = 0 spherical 
Bessel function of the first kind jo{z) = (sin z)/z. Defining the Fourier-Bessel transforms of 
momentum-space Salpeter component p^ip) and kinetic-energy contribution E{p) 922 (p) by 

dpp‘^ jo{pr)p 2 {p) , T{r) = dpp‘^ jo{pr) E{p) ip 2 {p) , (9) 

0 0 

Eq. (8) — the Bethe-Salpeter quintessence relevant here — becomes in conhguration space 

T(r) -I- V (r) ip{r) = 0 . 

Thus, anticipating that (p{r) has no zeros, for M = 0 the central potential in question reads 

V(r) = . (10) 

(p{r) 

N.B.: From Eqs. (9) and (10), both T(r) and V (r) approach, for large m, their trivial limits 
E{p) -)■ m T{r) - > mpir) V{r) - > —m . 

m—>-oo m^oo m—yoo 

4 Salpeter Amplitudes of Light Pseudoscalar Mesons 

At this stage, the only ingredient to the application of our inversion approach still lacking is 
the Salpeter amplitude of the bound states in the focus of our interest. In order to distill, by 
a heuristic line of argument, at least some rough idea of the shape of the Salpeter amplitude 
describing the light pseudoscalar meson, we take advantage of a relationship, proven within 
the context of Dyson-Schwinger equations [16,17], between this Salpeter amplitude, on the 
one hand, and the quark mass function in the dressed quark propagator, on the other hand. 
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4.1 Preliminary Definitions 

For simplicity of notation, let us introduce the bound-state vertex function r(p, P), derived 
from the Bethe-Salpeter amplitude <F(p, P) by removal of both fermion propagators Si{pi): 

T{p, P) = S^\p,) $(p, P) S^\-p2) ^ Hp, P) = ^i(pi) T{p, P) S2{-P2) . 

Using this quantity where convenient, the homogeneous Bethe-Salpeter equation (1) reads 

i 


r(p,p) = 


(27r) 


d^qK{p,q,P)^q,P) . 


Any exact fermion propagator S'(p) arises as a solution of the fermion Dyson-Schwinger 
(or gap) equation and, by Lorentz covariance and parity conservation, has to be of the form 

i 


S{p) = 




e 1 0 , 


AijP')— Bijp) + \e — MijP')+ ie ’’ 

where A{p^) and B{p^) are two real Lorentz-scalar functions which may be reinterpreted as 
the mass and wave-function renormalization functions of the respective fermion by dehning 


M(p^) = 


P(p2 


^Z{p^)B{p^) , Z{p^)^ 


AijA 


4.2 Results from Euclidean-Space Dyson—Schwinger Formalism 

For sound reasons [18], the apparatus of Dyson-Schwinger equations is typically developed 
in Euclidean space, with metric = 5^,^. Hence, for clarity, we identify in the following all 
Euclidean-space variables by underlining. Ignoring overall quark-flavour factors, the vertex 
function T{k, P) for a generic pseudoscalar meson P, with leptonic decay constant /p, reads 

r(fc,P) = 75 [iP(i,£) + ■■■], 

where E{k, P) is the dominant Dirac component, and the dots indicate the contributions of 
subleading [17] Dirac components. We can achieve our goal by combining two observations: 

1. In the chiral limit, the renormalized axial-vector Ward-Takahashi identity relates the 
Dirac component E{k, P) for P = 0 and the quark self-energy function P(£) [16,17]: 

/pP(£0) = P(fc2) . 

In the chiral limit, the vertex function of a (because of its Goldstone nature) massless 
(P^ = 0 ) pseudoscalar meson thus reads, in the center-of-momentum frame (P = 0 ), 


T{k,0) = 75 [iE{k,0) H-] = 75 


L/p 


Reinstalling both quark propagators yields the associated Bethe-Salpeter amplitude 
<!>(*:,0) = S(k) m, 0 ) S(k) cc 75 + • ■ ■ . 

2. We deduce the form of ^{k, 0) from an explicit solution [17] for the quark propagator: 

• The wave-function renormalization is usually very close to unity, i.e., Z{1^) ^ 1. 
Letting, for simplicity, Z{I^) ^ 1, i.e., M(fc^) ^ P(fc^), yields our starting point 


<h(fc, 0 ) oc 




■ 75 H- 


( 11 ) 
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• The Dyson-Schwinger equations constitute an infinite tower of coupled integral 
equations for the infinite number of n-point Green functions of a given quantum 
field theory. This system of equations relates every n-point Green function to, at 
least, one n'-point Green function with n' > n. Accordingly, finding solutions to 
such infinite hierarchy of Dyson-Schwinger equations requires, in the first place, 
the formulation of a manageable problem by truncation to a small finite number 
of relations for the low-order Green functions. Each higher-order Green function 
demanded as input by the truncated system of relations can only be modelled in 
accordance with all its general features expected from the quantum field theory. 
Glearly, both internal consistency and physical meaningfulness require any such 
truncation to be compatible with, at least, each of the identities in the totality of 
Ward-Takahashi identities encoding the symmetries of the underlying quantum 
field theory that proves to be indispensable for the problem under consideration 
[19]. A truncation scheme claimed to preserve the axial-vector Ward-Takahashi 
identity is the rainbow-ladder truncation, requiring a couple of approximations: 

— The exact (dressed) quark-gluon vertex function is replaced by its tree-level 
approximation; this change is usually dubbed as “rainbow approximation.” 

— At each occurrence, the Bethe-Salpeter kernel K is reduced to (an iteration 
of) its lowest-order perturbative contribution, single-gluon exchange, which 
results in the “ladder approximation” of any Bethe-Salpeter-type equation. 

— The exact (dressed) gluon propagator is replaced by its free approximation. 

— An appropriate effective coupling replacing the product of strong couplings 
in the ladder kernel thus obtained takes care of all phenomenological issues. 

A “renormalization-group-improved” variant [17] in the class of rainbow-ladder 
truncation models is specified by two crucial properties of the effective coupling: 

— In the infrared, —)■ 0, it shows the pronounced enhancement suggested by 

solutions of the Dyson-Schwinger equation for the exact gluon propagator. 

— In the ultraviolet, —)■ cxd, it approaches the perturbative behaviour of the 

strong fine-structure coupling and reproduces trivially asymptotic freedom. 


In the chiral limit, obtaining as solution of the Dyson-Schwinger equation for the full 
quark propagator nonvanishing dynamical quark masses M{k^) is a consequence and 
hence a signal of dynamical breakdown of chiral symmetry. In that model of Ref. [17], 
the one-loop behaviour of the quark mass function M{k^) for large Euclidean relative 
momenta k may be cast into a shape that involves the anomalous mass dimension 7 ^: 


lim M{k^) oc 

>-oo 


1 


7m 


12 

33-2Nf ■ 


Putting things together, the Bethe-Salpeter amplitude (11) behaves, in the ultraviolet, like 


«>(i. 0 ) oc 


75 • 


Ignoring the logarithmic correction, the Bethe-Salpeter amplitude is thus characterized by 


<h( 0 , 0 ) oc 


M(0) 


75 


lim ^{k, 0 ) oc -4 75 . 


>-oo 


k^ 


Mimicking the definition of the Salpeter amplitude by integration w.r.t. Euclidean time 
we find that Salpeter amplitudes of light pseudoscalar mesons fall off like |fc| for large \k\. 
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4.3 Implementation in Minkowski Space 

In the course of the above sketched expedition or raid into the jungle of Euclidean space, we 
could capture, as our main loot, two hints that might be of help in our quest for a justihable 
ansatz for the potential behaviour of the Salpeter amplitudes of light pseudoscalar mesons. 

• For large p, this amplitude 4>{p) decays proportional to the inverse third power of |p|; 

OC —rr for p —)■ CXD . 

|Pr 



• At the origin of three-momentum space, the amplitude 4>{p) can be taken to be hnite: 

1992(0)1 < 00 . 


Any justihable p dependence of the radial factor (p 2 ip) of the Salpeter component ip 2 {p) 
should conform to these boundary conditions. Hence, upon introducing a parameter p with 
the dimension of mass, our hope for analytic manageability prompts us to adopt the ansatz 

_ _ CXD 

^2(P) = 4 y ^ + p2)3/2 ’ h>0, \\^2f = j dpp^\^2{p)\'^ = l ■ (12) 

0 


The Fourier-Bessel transform (p{r) of such ansatz for the component 992 (p) of sole relevance 
for the amplitude 0(p) is nothing but the modihed Bessel function Kn{z) [8] of order n = 0: 


(p(r) 


4v^ 

71 


Koipr) , 


/i > 0 , 



(13) 


The behaviour of this ansatz in both momentum and conhguration space is shown in Fig. 1. 
Whereas 992(p) behaves smoothly, 99(r) diverges logarithmically at r = 0: 992 ( 0 ) = 4/ 


<^2(p) -1 0 , 

p—yoo 


(p(r) 


4v^ 




71 


ln(/ir) 


-)■ CXD 


r—>-0 


99(r) -)■ 0 . 

r—>-OD 


5 Configuration-Space Radial Potential by Inversion 

Our not too ambitious aim is to perform the inversion of the Bethe-Salpeter problem posed 
by the tentative Bethe-Salpeter amplitudes (12) or (13) as far as possible along an analytic 
path. That is to say, for as many different values of the free quantity m/p as manageable we 
intend to derive the analytical relationship between this ansatz for the pseudoscalar-meson 
solutions and the underlying Bethe-Salpeter kernel and to extract the interaction potential 
V (r) as a closed-form expression. It is straightforward to identify (a few) choices of m/p for 
which this task is easily accomplishable. For instance, for m = 0 (Subsec. 5.1) or for m = p 
(Subsec. 5.2.1), in which case the kinetic term E{p) (p 2 {p) entering V (r) via T(r) reduces to 

E{p)<p2{p)oc ^ (14) 

since the denominator of the p dependence of ^p 2 ip) equals [E{p)]^. In general, however, one 
has to content oneself with a numerical construction of the potentials V (r) (Subsec. 5.2.2).^ 


^ A rather condensed preliminary account of some of the results in this section may be found in Ref. [20]. 
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(a) (b) 

Figure 1 : Salpeter amplitude (2) for pseudoscalar mesons: independent component relevant 
for, at least, one (Fierz-invariant) Lorentz structure F(8)F = ^ (7/^(8)7^+75075 — 1 < 8 ) 1 ) of the 
interaction kernel in (a) momentum space and (b) configuration space, given in units of the 
parameter p, i.e., for p = 1 . In momentum-space representation, (^2(p) = 4 
approaches for p —)■ 0 the hnite value <^2(0) = 4/y^ 2.256785 .... In conhguration-space 

representation, on the other hand, <p(r) = (4 \/2/n) Ko{r) necessarily encounters for r —)■ 0 
the logarithmic singularity of the modihed Bessel function Ko{z): ip{r) —)■ —(4 -s/^/ti) ln(r). 


5.1 Bound-state constituents of vanishing mass: massless quarks 


In the ultrarelativistic limit, things become very simple. For m = 0 and therefore E{p) = p, 
the Fourier-Bessel transform T(r) of the kinetic term involves the modihed Bessel function 
In{z) [8] (of order n = 0,1) and the modihed Struve function L„(;2) [8] (of order n = —1, 0): 

r(r) = ^ [/o(pr) + pr/i(pr) -/irL_i(/ir) -Lo(/ir)] . 

r 

So, a potential (10) yielding massless pseudoscalar bound states of massless constituents is 

^ ^ ^ /xrL_i(/rr) + Lo(/rr) - Ipjiar) - pr/i(/ir) 

2 r Ko{pr) 

As depicted in Fig. 2, this potential is characterized by a logarithmically softened Coulomb 
singularity at the origin r = 0 and a confining rise beyond bounds for large distances r, i.e.. 


V{r) 


TT 


2rln(pr) r->-o 


-00 


V{r) 


TT r’^ 


exp(/ir) 


-)■ 00 


The negative portion of V(r) for small distances r, in cooperation with the strongly peaked 
r-dependence of (p(r) near r = 0, serves to counterbalance the positive contributions to the 
bound-state mass M arising from the kinetic term and from the positive portion of V(r). It 
is therefore crucial for maintaining the desired masslessness of the bound state, i.e., M = 0. 
In Subsec. 5.2.2, the above r dependence will prove to form a prototype for the behaviour of 
the potential V(r) for any mass m of the bound-state constituents in the region 0 < m < p. 
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Figure 2: Configuration-space potential V (r) deduced by inversion of the Salpeter equation 
(3) with the Lorentz structure F^F = ^ (7;i(8)7^+75(8)75 — 1 ( 8 ) 1 ) of the interaction kernel for 
the ansatz </? 2 (p) oc (p^-|-l)“^/^ of the relevant component of the Salpeter amplitude 0(p) in 
momentum space that is expected to describe massless (M = 0) pseudoscalar bound states 
of zero-mass (mi ,2 = 0) constituents: V{r) = vr [r L_i(r)-|-Lo(r)—/o(r)—r /i(r)]/[2 r K^ij)] 
exhibits a singularity (r In r)“^ at r = 0 and a confining rise to infinity for large distances r. 


5.2 Bound-state constituents of non-zero mass: massive quarks 

5.2.1 Case 0 < m = fi 

If the non-vanishing common mass m of the two bound-state constituents is precisely equal 
to the free “smoothing” parameter p, that is, for 0 < m = p, the kinetic term E{p) (p 2 {p) is, 
from Eq. (14), proportional to the Fourier transform of the Yukawa shape exp{—z)/z] thus, 
its Fourier-Bessel transform T(r) is inevitably of Yukawa form, with m as slope parameter: 


T{r) 


2\/2n? 


exp(—mr) . 


r 
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Accordingly, under the circumstances discussed above, a potential that leads, as solution of 
the Salpeter equation, to massless pseudoscalar bound states of massive constituents reads 

^ ^ TT exp(-mr) ^ 

2 rKo{mr) 

Similarly to the case of zero-mass constituents inspected in Subsec. 5.1, the above potential 
has a logarithmically softened Coulomb singularity at the origin r = 0, where the Coulomb 
singularity reflects the Yukawa-type shape of T(r), whereas the logarithmic softening again 
derives from the (singular) behaviour of the modihed Bessel function Kq{z —)■ 0) ~ — ln(; 2 ); 
in contrast, for large inter-constituent separation r the potential approaches zero like 

^ T ^ ^ 10 . 

r^O 2rln(?7ir) r^o r-^-oo y 2r r—>-oo 

Such singular V (r) shape (Fig. 3) cancels exactly the contribution of the kinetic term T(r); 
this kind of nonconfining behaviour will turn out to be generic for arbitrary ratio m/fa > 1. 

5.2.2 Case 0 < m 

For arbitrary non-vanishing values of the common mass m of the bound-state constituents, 
i.e., for 0 < m 7 ^ /i, expecting the resulting potential to be expressible in closed form would 
betray an incommensurately high degree of optimism. In general, we have to be satished by 
extracting the potential V (r) by numerical evaluation of Eq. (10). In order to get some idea 
about the emerging regularities. Fig. 4 depicts the outcomes of this extraction for a number 
of selected representative values of m. Inspection of these Endings leads us to conclude that 
the crucial quantity determining the behaviour of V (r), at least for large separation r of the 
constituents, is the relative magnitude of their mass m and the wave-function parameter /r: 

• At the origin r = 0, the potential V (r) seemingly displays, irrespective of the value of 
the constituents’ mass m, a logarithmically softened Coulomb singularity of the form 

V{r) - > ——-- > —oo . 

r^-o 2rln(/ir) r^o 

• For masses m smaller than the parameter fa, 0 < m < fi, the potential V (r) increases, 
for large r, without limits, and constitutes thus a manifest realization of confinement. 

• For masses m not less than the parameter fa, 0 < fi < m, the potential V (r) tends, for 
rising r, towards a nonpositive constant that, in the limit of large m, approaches —m: 

lim V (r) - > —m . 

r^oo m^oo 


6 Summary of Findings 

Combined Dyson-Schwinger-Bethe-Salpeter studies show that the bound-state amplitude 
of a light pseudoscalar meson decreases, for very large relative distances between the bound 
quark and antiquark, approximately like a power law. In the present analysis, we addressed 
the seemingly innocent question how interaction potentials that — when being fed into the 
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Figure 3: Configuration-space potential V (r) deduced by inversion of the Salpeter equation 
(3) with the Lorentz structure F^F = ^ ( 7 ;i( 8 ) 7 ^+ 75 ( 8)75 —1(8)1) of the interaction kernel for 
the ansatz </? 2 (p) oc (p^-|-l)“^/^ of the relevant component of the Salpeter amplitude 0(p) in 
momentum space that is expected to describe massless (M = 0 ) pseudoscalar bound states 
of constituents with non-vanishing masses mi ^2 = 1: V{r) = —7rexp(—r)/[2r iFo(?")] shows 
for r —)■ 0 the same logarithmically softened Coulomb singularity (r lnr)“^ as the potential 
found for zero-mass constituents (cf. Fig. 2) but, in distinct contrast to the behaviour of the 
latter, approaches for r —)■ cxd a finite value, 0 , and is, accordingly, a nonconfining potential. 


instantaneous Bethe-Salpeter formalism — yield bound states with asymptotic power-law 
decrease might look like. The resulting potential exhibits presumably unexpected features. 
For increasing interquark separation, the potential is monotone rising. It starts at a sharply 
peaked Coulomb-like singularity at the origin but flattens off at intermediate distances. For 
large distances, it rises, for both quarks sufficiently light, rather steeply to infinity, whereas, 
for all heavier quarks, it approaches a nonpositive hnite value; the borderline for the change 
of behaviour is drawn by the size of a mere model parameter which, however, in a genuinely 
fundamental treatment will result from the basic parameters of quantum chromodynamics. 



12 


0 2 4 6 8 10 12 14 16 



16 

12 

8 

4 

0 

-4 

-8 


-12 


-16 


Figure 4: Configuration-space potential C(r) deduced by inversion of the Salpeter equation 
(3) with the Lorentz structure F^F = ^ ( 7 ;i( 8 ) 7 ^+ 75 < 8)75 —1<8)1) of the interaction kernel for 
the ansatz </? 2 (p) oc of the relevant component of the Salpeter amplitude in 

momentum space that is expected to describe massless (M = 0) pseudoscalar bound states 
of constituents with arbitrary masses 0 < mi = m 2 = m. In two exceptional cases, a closed 
expression for V (r) may be obtained: V(r) = tt [r L_i(r)-|-Lo(r)—/o(r)—r /i(r)]/[2 r -K’o(r)] 
for m = 0 (dotted line, see Fig. 2) and V (r) = —tt exp(—r)/[2 r iFo(r)] for m = 1 (solid line, 
see Fig. 3); for other masses m exemplifying the general case, m = 0.35 (short-dashed line), 
m = 0.5 (long-dashed line), and m = 2 (dot-dashed line), V (r) has been found numerically. 
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